Joule heating generated by spin current through Josephson junctions 



S. Takahashi,* T. Yamashita, T. Koyama, and S. Maekawa 
Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 



o 
o 
o 

(N 

■4— » 

o 
O 



o 
o 

Of 



-a 
c 

o 
o 



> 
o 

CO 
CN 

o 

o 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



H. Imamura 

Graduate School of Information Sciences, Tohoku University, Sendai 980-8579, Japan 

We theoretically study the spin-polarized current flowing through a Josephson junction (JJ) in 
a spin injection device. When the spin-polarized current is injected from a ferromagnet (FM) in 
a superconductor (SC), the charge current is carried by the superconducting condensate (Cooper 
pairs), while the spin-up and spin-down currents flow in the equal magnitude but in the opposite 
direction in SC, because of no quasiparticle charge current in SC. This indicates that not only the 
Josephson current but also the spin current flow across JJ at zero bias voltage, thereby generating 
Joule heating by the spin current. The result provides a new method for detecting the spin current 
by measuring Joule heating at JJ. 



Spin transport through a nonmagnetic metal has at- 
tracted much interest in magnetic nanostructures. In 
the tunnel junctions consisting of a ferromagnet (FM) 
and a normal metal (N) or superconductor (SC), the 
tunnel current driven from FM is spin-polarized |jj and 
creates a nonequilibrium spin population in N or SC 
P,PJ- Recently, there has been a number of experiment 
on suppression of superconductivity by injection of spin- 
polarized electrons using tunnel junctions of a high-T c 
SC and a ferromagnetic manganite @,[|,U ■ 

A double tunnel junction with a thin layer of SC sand- 
wiched between two FM electrodes is a unique system 
to investigate a novel magnetoresistive effect caused by 
nonequilibrium spin population in SC. When the thick- 
ness of SC layer is much smaller than the spin diffusion 
length, the spin population depends strongly on the rela- 
tive orientation of magnetizations in FMs. In the antipar- 
allel alignment of magnetizations, the spin population is 
accumulated in SC and reduces the superconducting gap 
A with increase of tunneling currents, while in the par- 
allel alignment there is no such effect because of the ab- 
sence of spin accumulation in SC j?J . In contrast, the spin 
current flows in SC in the parallel alignment, but not in 
the antiparallel alignment ||] . Although several methods 
for detecting the spin accumulation has been presented, 
there is little for the spin current. In this paper, we pro- 
pose a new method for detecting the spin current flowing 
in SC using a Josephson junction (JJ). We show that 
the spin current flows through JJ even at zero bias volt- 
age (Vj = 0), thereby generating Joule heating at JJ. 
Since there is no quasiparticle charge current across JJ 
at Vj =0, the Joule heating generated at Vj = is a 
direct signature of the spin current flowing in SC. 

We consider a spin injection device made of FM and 
SC separated by thin insulating layers as shown in Fig. 1. 
The central part of the junction forms the Josephson 
junction, which is sandwiched between two FMs. The 
left and right SC (FM) are made of the same SC (FM). 
The magnetization of the left FM points up and that of 
the right FM is either up or down. The applied bias 



current ij n j flows through the junctions of resistances Rj 
and Rt with the voltage drops Vj and Vr, respectively. 

We calculate the tunneling current using a phenomeno- 
logical tunneling Hamiltonian that describes the transfer 
of electrons from one electrode to the other. If SC is in 
the superconducting state, it is convenient to rewrite the 
electron operators ako- in SC in terms of quasiparticle 
operators 7k CT appropriate to the superconducting states, 
using the Bogoliubov transformation 
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where |u k | 2 = |(l + £ k /£ k ), 



is an operator which annihilates a Cooper pair, while S' 
creates one, and = + A 2 ] 1 / 2 is the quasiparticle 
dispersion of SC, £ k being the one-electron energy relative 
to the chemical potential of the condensate and A being 
the isotropic superconducting gap. 

From the Fermi's golden rule result, the spin- 
dependent tunnel currents Ii„ across the zth junction are 
expressed as [Q 
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Here, Gi„ (i — 1,2) is the tunnel conductance for elec- 
trons with spin a when SC is in the normal state. The 
quantity J\f is given by the usual expression 



AA(Vt) 



V s (E)[f (E-eV T )-f (E)]dE, (2) 



where T> S (E) = Re [\E\/\/E 2 - A 2 J is the normalized 
BCS density of states and fo(E) the Fermi function. The 
quantity <Sj is the quasiparticle spin density accumulated 
inSC 
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FIG. 1. Spin injection device consisting of a Josephson 
junction (SC-I-SC) which is sandwiched between two ferro- 
magnets (FM). The arrows indicate the injection current, pair 
current, and quasiparticle current in SCs, and the Josephson 
current (solid arrow) and spin current (dotted arrow) across 
JJ, in the case of parallel alignment of magnetizations in FMs 
whose spin polarization is 100 %. 



terms describe the phase coherent (Cooper pair) tunnel- 
ing. The usual Josephson effect is associated with the 
sin ip term. Using the golden rule formula, we have the 
spin-dependent tunnel current of quasiparticles 



x [f 1(7 (E)-f 2a (E + eVj)]dE, 



(9) 



with which the quasiparticle charge current 7 qp and the 
spin current J sp i n across JJ are written as 



tqp 
*spin 



/smn = 4(Vj) 



(10) 
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The phase coherent tunneling terms are obtained as 
1 V S (E)V S (E + eV 3 ) 



J - A2 



dE- 



\E(E + eVj)\ 
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(12) 



1>s(E)[f ir (E) - fn(E)]dE, 



(3) 



where fi a (Ek) — (7^. (T 7kcr)j is the distribution function of 
quasiparticles with energy E^ and spin a in the ith SC. 
In Eq. (1), we neglected the charge imbalance (Q*) since 
it has no relevance to the spin-dependent effect 

In the limit of vanishing spin- flip scattering in SC, the 
spin-up and spin-down currents are treated as indepen- 
dent channels. The conservation of the total charge cur- 



rent 



/ and the spin current / s . 



across the left and right junctions yields the relations 
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[{P 1 + P 2 )Af-{S 1 -S 2 )}, 



(4) 
(5) 



where 1/R T = G iT + G a and P t = (G iT - G U )/(G^ + 
Gii). Pi represents the degree of spin-polarization of FM. 
In the parallel alignment, P\ = P 2 = P and S\ = —52 = 
5, so that 



/inj " ei? T 



(TV - PS) 



/spin ~ eR T 



(PAf-S). (6) 



In the antiparallel alignment, Pi 
S 2 = PAf, so that 



-P 2 = P and Si 



J inj = (eRr)- 1 (1 - P 2 ) TV, J spin = 0. 
The tunnel current across JJ has the form 
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^qp(Vj) + iji(Vj) sin<^ + Ij 2 (Vj) cos(^, 



(7) 



(8) 



for a bias voltage Vj and a phase difference ip of the 
gap parameters on the two side of J J. The first term de- 
scribes the quasiparticle tunneling, and second and third 
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(13) 



When the thickness of SC is much smaller than the spin 
diffusion length, the distribution of quasiparticles is spa- 
tially uniform in SC. Then, the distribution function fi a 
is described by /o, but the chemical potentials of the spin- 
up and spin-down quasiparticles are shifted oppositely by 
from the equilibrium one to create the spin density; 

fa(E) = f (E - 6m), fn(E) = f (E + 6m). (14) 

The gap parameter A in SCs is determined by fi a 
through the BCS gap equation 
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dE = 0, 



(15) 



where Ao is the gap at T = in the equilibrium state 

(Sm = q). 

In the following we calculate the tunnel current at zero 
bias voltage (Vj =0). In this DC Josephson cur- 

rent flows across JJ unless the bias current exceeds the 
Josephson critical current J c = lj(0). From Eq. (|Tl|), we 
notice that /spin in the parallel alignment becomes finite 
even at Vj = if S^i and bp, 2 take nonzero values of 
different sign = 5 Hi = —6fj, 2 ), while I sp i n in the an- 
tiparallel alignment is zero because of Sfj, = Shi = &\x 2 . 
The currents IJ2 and I qp vanish at Vj = irrespective of 
the value of 6 Hi. These facts indicate that the Josephson 
current as well as the spin current flow across JJ at Vj = 
in the parallel alignment, while only the Josephson cur- 
rent flows in the antiparallel alignment. 
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FIG. 2. Joule heating W as a function of temperature at 
zero bias (Vj = 0) in the parallel alignment of magnetizations 
in FMs. Inset shows half the spin splitting of the chemical 
potential of quasiparticles in SC. 



We consider, for simplicity, the case where FM is a 
half-metal with 100 % spin polarization (P = 1), in which 
the spin current is equal to the injection current, /i n j = 
-fspin, because only the spin-up electrons are injected into 
SC. Figure 1 shows how the current 7 in j = jf is injected 
from the left FM in SC and how the spin and charge 
currents flows in SCs and pass through JJ in the parallel 
alignment. In SC, the charge current is carried by the 
superconducting condensate (Cooper pairs) alone; there 
is no charge current carried by quasiparticles (I^ p + I^ p — 
0) , so that the spin- up and spin-down currents flow in the 
opposite direction with the value of J? = —Iq P = 1^/2. 
In other words, the half of /T is injected to be the spin- 
up current 1^ in SC, and the other half of I' injected 
in SC is combined with the backfiow spin-down current 
J^p to form the pair current I pa ir- It is remarkable that 
the spin current is able to pass through JJ at zero bias 
(Vj = 0). This is because the voltage drop across J J 
is positive for the spin- up current, while it is negative 
for the spin-down current, resulting in no net voltage 
drop across JJ. Therefore, the spin current as well as the 
Josephson current flow through the JJ at Vj = 0, the 
magnitude of which is equal to ij n j in the half-metallic 
case. 

The most striking prediction of the present theory 
is that the spin current across JJ is accompanied by 
Joule heating W at zero bias voltage (Vj = 0). The 
Joule heating W generated at Vj = is given by 
W = /spin (<fyti — Sfj.2) /e. In the parallel alignment 



(5(i = 5 hi = -5(j, 2 ), 



W = 2I spin 5fi/e, 



(16) 



while in the antiparallel alignment W = because of 
^spin = and Sfi = 5fj,2 = 8fj,i, indicating that the 
Joule heating generated by spin current has strong spin- 
dependence. 



Figure 2 shows the Joule heating W normalized to 
Wo = Ap/e 2 i?j as functions of temperatures for the 
injection current I-mj/JcO = 0.1 and 0.2 in the parallel 
alignment |10|, where J c o = 7rAo/(2ei?j) is the Joseph- 
son critical current. As the temperature is lowered, W 
increases monotonically. This is due to the increase of 
the spin-splitting of S/j, (see the inset) in the current bias 
mode. A slight decrease of the critical temperature and 
a depression of W near T c in the curve of Ii n j/Jc0 = 0.2 is 
caused by the pair breaking effect by spin accumulation 
in SC. The power of Joule heating W generated by spin 
current is estimated by evaluating the value of Wo, since 
W is comparable to Wo- If one uses the values of an area 
resistance RjA = 100 Sl^m 2 and A = 0.3 meV (Al), we 
have Wo /A — 90 mW/cm 2 per area of J J. Therefore, W 
is fairly large for observing the Joule heating experimen- 
tally. 

In conclusion, we propose a new method for detecting 
the spin current by measuring the Joule heating gener- 
ated at JJ. This is attributed to the fact that the spin- 
up and spin-down currents flow in the opposite direction 
in SC. For measuring the Joule heating, it may be im- 
portant for the tunnel junctions to satisfy the condition 
Rj ^> in order to make the Joule heating W domi- 
nant compared with the injection power (I? n -Rrr). 

We are grateful to M. Johnson and J. S. Moodera 
for valuable discussions. This work is supported by a 
Grant-in-Aid for Scientific Research Priority Area for 
Ministry of Education, Science and Culture of Japan, 
NEDO Japan, and by the supercomputing facilities in 
IMR, Tohoku University. 



[1] R. Meservey and P. M. Tedrow, Phys. Rep. 238, 173 
(1994). 

[2] M. Johnson and R. H. Silsbee, Phys. Rev. Lett. 55, 1790 

(1985); Phys. Rev. B 37, 5326 (1988). 
[3] M. Johnson, Appl. Phys. Lett. 65, 1460 (1994). 
[4] V. A. Vas'ko et al., Phys. Rev. Lett. 78, 1134 (1997). 
[5] Z. W. Dong et al, Appl. Phys. Lett. 71, 1718 (1997). 
[6] N.-C. Yeh et al, Phys. Rev. B 60, 10522 (1999). 
[7] S. Takahashi, H. Imamura, and S. Maekawa, Phys. Rev. 

Lett. 82, 3911 (1999); J. Appl. Phys. 85, 5227 (2000). 
[8] S. Takahashi et al. Physica C, (2000). 

[9] B. D. Josephson, Phys. Lett. 1, 251 (1962); M. Tin- 
kham, Introduction to Superconductivity (McGraw-Hill, 
New York, 1996). 
[10] To avoid a spurious logarithmic divergence arising from 
convolving the BCS density of states with an isotropic 
gap in Eq. ([n]), we introduce a small broadening in the 
BCS density of states using Dynes's formula T>s{E) = 
Re[(_B — 17) / (E — ij) 2 — A 2 ] with broadening parame- 
ter 7 much less than A. 



3 



